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Abstract 

In this review we demonstrate how the algebraic Bethe ansatz is used for the 
calculation of the energy spectra and form factors (operator matrix elements in the 
basis of Hamiltonian eigenstates) in exactly solvable quantum systems. As examples 
we apply the theory to several models of current interest in the study of Bose- 
Einstein condensates, which have been successfully created using ultracold dilute 
atomic gases. The first model we introduce describes Josephson tunneling between 
two coupled Bose-Einstein condensates. It can be used not only for the study 
of tunneling between condensates of atomic gases, but for solid state Josephson 
junctions and coupled Cooper pair boxes. The theory is also applicable to models of 
atomic-molecular Bose-Einstein condensates, with two examples given and analysed. 
Additionally, these same two models are relevant to studies in quantum optics. 
Finally, we discuss the model of Bardeen, Cooper and Schrieffer in this framework, 
which is appropriate for systems of ultracold fermionic atomic gases, as well as being 
applicable for the description of superconducting correlations in metallic grains with 
nanoscale dimensions. In applying all of the above models to physical situations, 
the need for an exact analysis of small scale systems is established due to large 
quantum fluctuations which render mean-field approaches inaccurate. 



*email: jrl@maths.uq.edu.au 
^email: hqz@maths.uq.edu.au 



1 



Contents 

ll Introducti 



4 



ion 



2 Model Hamiltoniansl 6 



2.1 A model for two Josephson coupled Bose-Einstein condensates! 7 



2.2 A model for homo-atomic-molecular Bose-Einstein condensates! 7 



2.3 A model for hetero-atomic-molecular Bose-Einstein condensates! 

2.4 The reduced BCS mo 



3 Quantum inverse scattering methodl 9 



3.1 The quasi-classical limitl 13 



3.2 Examples of realisation of the Yang-Baxter algebra! . 14 



3.3 Z-graded structure and Z-graded realisations of the Yang-Baxter algebra! . 15 

3.4 Examples of Z-graded realisations! * . . . . 17 



Algebraic Bethe ansatz method of solution! 18 



4.1 Extension to Z-graded realisations 21 



5 Scalar products of states! 23 



5.1 The quasi-classical limitl 26 



6 Exact solution of the models] 27 



6.1 Solution for the model of two Josephson coupled Bose-Einstein condensates! 27 



6.2 Solution for the model of homo-atomic-molccular Bosc-Eiiistcin condensates! 28 



6.3 Solution for the model of hetero-atomic-molecular Bose-Einstein condensates! 29 



<L4 Solution for the reduced BCS modell 31 



Exact calculation of form factors! 34 



7. 1 Form factors for the model of two Josephson coupled Bose-Einstein condensates! 34 

7.2 Form factors for the models of atomic-molecular Bose-Einstein condensates! 3£ 



7.3 Form factors for the reduced BCS modell 42 



8 Conclusion 44 



2 



1 Introduction 



Exactly solvable models of quantum mechanical systems provide an important insight 
into the nature of quantum physics, with the simple harmonic oscillator and the non- 
relativistic hydrogen atom serving as the archetypal examples. One method to solve 
these models is to exploit an underlying algebraic structure, well known to be the Lie 
algebra gl(3) for the harmonic oscillator and so(4) for the hydrogen model PQ. In this 
approach, the Lie algebraic structure plays the role of generating states of the system 
while at the same time providing state labels (or quantum numbers). A celebrated exact 
solution of a quantum many-body model is that for the one- dimensional Heisenberg (spin 
1/2) chain, due to Bethe |2j. Out of this work grew the concept of the Bethe ansatz for 
the construction of the eigenvectors for an exactly solvable Hamiltonian. In adopting this 
method, a general possible form for an eigenvector is assumed, dependent on several free 
parameters. Constraints are then determined for the parameters which ensure that this 
vector is an eigenvector of the Hamiltonian. The constraint equations are referred to as 
the Bethe ansatz equations of the model. 

Motivated by Bethe's work the field of exactly solvable models flourished during the 
1960s led by McGuire [3], Lieb [I], Sutherland [5], Yang jHj and Baxter [7j, amongst many 
others. Out of this activity arose the Yang-Baxter equation, the solution of which pro- 
vides a sufficiency condition to construct a model which is exactly solvable (applicable 
to one- dimensional quantum spin chains, including quantum field theories as the lattice 
spacing goes to zero, and classical two-dimensional lattice systems) jSJIH]. A fundamental 
feature of the Yang-Baxter equation is that it can always be used to construct a family 
of mutually commuting matrices, known as transfer matrices, which facilitates the appli- 
cation of the Bethe ansatz. The method of the Bethe ansatz can take a variety of forms, 
commonly known as the co-ordinate, analytic, functional and algebraic forms. It is this 
latter approach that will be the focus of our work here, as this is the most appropriate to 
serve our requirements. 

The algebraic formulation of the Bethe ansatz, and the associated quantum inverse 
scattering method, was primarily developed by the group of mathematical physicists in 
St. Petersburg (TU1 HU C21 ECU El- Its applicability extends beyond the study of one- 
dimensional spin chains, quantum field theory and two-dimensional lattice models to 
systems of correlated electrons [H^, conformal field theory [Ug^, as well as precipitating the 
notion of quantum algebras (deformations of universal enveloping algebras of Lie algebras) 
|17| H%1 IT9"1 l2*Uj . The main motivation for the algebraic formulation of the Bethe ansatz 
was not only for calculating the energy spectrum of a model, but to also accommodate 
the calculation of correlation functions. An initial step in this direction is to compute the 
form factors of an operator (not necessarily observable). Here the term "form factors of an 
operator" simply refers to the matrix elements of that operator in the basis of Hamiltonian 
eigenstates. Expectation values of observable operators and general correlation functions 
are expressible in terms of form factors through completeness relations. 

The study of correlation functions in the context of exactly solvable models has its 
origins in Baxter's corner transfer matrix method 7 j. Following on from this there has 
been a rich theory developed using ideas taken from affine quantum algebras, vertex oper- 
ators, integrable field theories, the off-shell Bethe ansatz and the Knizhnik-Zamolodchikov 
equation (2H 1221 1221 I2U 1213 1213 I2Z] , as well as the algebraic Bethe ansatz approach (e.g 
[E3 12H1 I2H1 EIB HH] ) • For the models we will study here the calculation of form factors will 
be undertaken through extensive use of the Slavnov formula [32] for the scalar products of 
Bethe eigenstates. The Slavnov formula provides an explicit determinant representation 
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for the scalar products. A refined proof of this result was given by Kitanine, Maillet 
and Terras [29J, using the notion of factorising the solutions of the Yang-Baxter equa- 
tion in terms of Drinfeld twists [33]. They applied this method to successfully compute 
form factors for the anisotropic (XXZ) Heisenberg chain and in a closely related 
work Korepin and Slavnov computed form factors for the quantum non-linear Schrodinger 
equation j30]. In both cases the results are valid for finite size systems, and thus this ap- 
proach is appropriate for applications to nanoscale systems. The results presented here 
are largely inspired by these works. 

The need to appeal to the exact solution of a model has been well illustrated in the con- 
text of the energy spectrum of metallic grains of nanoscale size. Experiments conducted 
by Ralph, Black and Tinkham (RBT) [3U EH] using single electron tunneling spectroscopy 
on aluminium grains with mean radii in the range 5-13 nm indicated significant parity 
effects due to the number of electrons in the system. The electron number remains fixed 
due to the large charging energy of the grains, which is a consequence of their small size. 
For grains with an odd number of electrons, the gap in the energy spectrum reduces with 
increasing size of the system, in contrast to the case of a grain with an even number of 
electrons, where a gap larger than the single electron energy levels persists. In the latter 
case the gap can be closed by a strongly applied magnetic field. The conclusion drawn 
from these results is that pairing interactions are prominent in these nanoscale systems. 
For a grain with an odd number of electrons there will always be at least one unpaired 
electron, so it is not necessary to break a Cooper pair in order to create an excited state. 
For a grain with an even number of electrons, all excited states have a least one broken 
Cooper pair, resulting in a gap in the spectrum. In the presence of a strongly applied 
magnetic field, it is energetically more favourable for a grain with an even number of 
electrons to have broken pairs, and hence in this case there are excitations which show no 
gap in the spectrum. 

A naive approach to describe these nanograms is to apply the theory of supercon- 
ductivity due to Bardeen, Cooper and Schrieffer (BCS) [36J. Indeed, the BCS model is 
appropriate for these systems but the associated mean-field treatment fails. There are two 
main reasons for this. First is because the BCS analysis makes use of the grand canonical 
ensemble whereas in the experiments the electron number is fixed. Second is because 
a mean-field theory approximates certain operators in the model by an average value. 
At the nanoscale level, the quantum fluctuations are sufficiently large enough that this 
approximation is invalid. In systems where the mean single particle energy level spacing, 
which is inversely proportional to the volume, is comparable to the bulk superconduct- 
ing gap, (as is for metallic nanograms,) it was thought that pairing interactions would 
not correlate any energy levels. This was conjectured by Anderson j37] on the basis of 
the BCS mean-field analysis, but the experiments of RBT show this to not be the case. 
Consequently, an exact solution was desired in order to clarify the issue. 

Remarkably, the exact solution of the reduced BCS model ("reduced" refers to the fact 
that only zero momentum Cooper pairs are considered and all couplings for scattering of 
Cooper pairs are equal) had been obtained and analysed many years earlier in a series of 
works by Richardson and Sherman [3BJ EH] using an approach equivalent to the co-ordinate 
Bethe ansatz. The motivation for their work was for application of pairing interactions in 
nuclear systems, which involve fixed particle number, and thus the BCS analysis referred 
above is not valid. However, the condensed matter physics community was unaware of 
this earlier work at the time the results of RBT were communicated. It was subsequently 
shown that theoretical results obtained through an analysis of the exact solution for the 
reduced BCS Hamiltonian were compatible with the experimental results of RBT [40J . 
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One of the most currently active fields is the study of Bose-Einstein condensates of 
ultracold atomic gases |HJ H2] • The Bose-Einstein condensed state is of a purely quan- 
tum mechanical nature and, in analogy with the phenomena of superconducting metallic 
nanograins discussed above, a mean-field analysis of small scale systems comprised of 
Bose-Einstein condensates is inadequate due to significant quantum fluctuations. While 
there are many novel physical properties to be explored in the study of Bose-Einstein 
condensates, there are three features that we will discuss here. The first is the phe- 
nomenon of Josephson tunneling between two coupled Bose-Einstein condensates. Recall 
that the Josephson effect was first proposed in relation to the tunneling of Cooper pairs 
through an insulating barrier separating two superconductors I44j . (A very infor- 
mative historical account is given in |4"5|.) It has been proposed as a means to couple 
qubits for the purpose of quantum computation jJHll^j. The experimental realisation of 
Bose-Einstein condensation in the atomic alkali gases provides a framework in which to 
observe macroscopic tunneling in a system with tunable couplings. An extensive account 
of this phenomenon can be found in [IS], which discusses in detail the canonical Joseph- 
son Hamiltonian (equivalent to a two site Bose-Hubbard model) for the description of 
this effect. It is not well known that this model is exactly solvable through the quantum 
inverse scattering method, which was established about a decade ago in the context of the 
discrete self-trapping dimer model [JOJ ED* • Below we show that a slightly more general 
model is also exactly solvable and we derive explicit exact form factors for the generalised 
model. 

The second aspect of Bose-Einstein condensation we will discuss is that of a condensate 
comprised of a coherent superposition of atomic and molecular states. This phenomenon 
has been predicted and studied by theorists (e.g., see j^H E21 EE1 EH ES]) and recently 
realised experimentally EH EH] • in particular, for the experiment of JSH] using 85 Rb 
atoms, which are converted into diatomic molecules via a Feshbach resonance, the sys- 
tem was prepared, allowed to evolve, and then a measurement made to determine the 
number of atoms in the system. By performing this procedure over different evolution 
times it was established that the expectation value for the number of atoms displayed an 
oscillatory behaviour, indicating that the state of the system was a quantum mechanical 
superposition of atomic and molecular states, as opposed to a classical mixture. The 
result is significant in that the state of the system is comprised of a superposition of two 
chemically distinct components. 

Finally, we will also analyse the reduced BCS model, which is relevant not only for 
metallic nanograins as described above, but also the study of ultracold fermionic atomic 
gases jSS]. As is well known, for an ultracold fermionic gas the Pauli principle prohibits 
all particles occupying the lowest energy level. The lowest possible energy of the system 
is obtained by filling the Fermi sea. However, in analogy with metals it is believed that 
fermionic gases should be able to form Cooper pairs, and as a consequence, undergo a 
phase transition at a suitably low temperature into a fermionic condensate [HD] • 

The aim of this exposition is to illustrate that the algebraic approach to the study of 
exactly solvable models is a rich and elegant theory with wide applicability. In particular 
we show how the theory applies to the systems of Bose-Einstein condensates and the 
reduced BCS model discussed above. Some of these results have already been communi- 
cated jHIl E21 EHI1 EH EH EH] , while other results we will present are new. In each case we 
will determine the energy spectrum, as well as the form factors for the computation of 
correlation functions, in terms of the Bethe ansatz solution. Certain correlation functions 
can in fact be deduced directly from the energy spectrum through use of the Hellmann- 
Feynman theorem [HZ1 EE]- Examples of this procedure applied to the models discussed 
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here can be found in |U3J EH EH] ■ Typically however, the form factor approach is required 
to build general formulae for expectation values and correlation functions. As a potential 
application of these results we point to the problem of quantifying entanglement in the 
theory of quantum information. The role of correlation functions in the characterisation 
of entanglement has been discussed in (EH EH EU [72] . 

Throughout, we have endeavoured to provide as much technical detail as possible for 
the benefit of non-experts. The exceptions are the Slavnov formula for the scalar product 
of states, the proof of which is beyond the scope of this review. For the proof we refer 
the interested reader to Also, the orthogonality of the Bethe eigenstates will not 

be proved. Details of this result can be found in [SO]- The format of the review is as 
follows. We begin in section 2 with a description of the four models we will examine. 
In section 3 we recall the basic features of the quantum inverse scattering method for 
the construction of exactly solvable models. While there already exist several excellent 
surveys of this approach [TU1 HU E21 EE El, we give a detailed account here in order to 
fix notations and conventions and make the review self-contained. The central aspect is 
the introduction of the Yang-Baxter algebra associated with the Lie algebra gl(2), which 
is a quadratic algebra. Several examples of realisations are given. We show that in a 
particular limit, called the quasi-classical limit, the Yang-Baxter algebra reduces to a Lie 
algebra, called the Gaudin algebra. Through a realisation of the Yang-Baxter algebra, 
the transfer matrix is constructed which leads to an exactly solvable model. We also 
discuss a natural Z-graded structure of the Yang-Baxter algebra which will be exploited 
in later constructions. Section 4 deals with the algebraic Bethe ansatz method in a general 
context for the determination of the spectrum of the transfer matrix. Section 5 presents 
the Slavnov formula for the scalar products of the states which arise in the algebraic 
Bethe ansatz method of solution. We also discuss how, through the use of the Slavnov 
formula, the form factors for the elements of the Yang-Baxter algebra can be obtained. 
Section 6 turns to calculating the explicit exact solutions for the models. Formulae for 
the energy spectrum are determined, which are parameterised in terms of the roots of the 
Bethe ansatz equations. Section 7 deals with the computation of form factors for each 
of the models introduced. In all cases it is necessary to first consider the solution to the 
inverse problem, which involves expressing a given operator in terms of the elements of 
the Yang-Baxter algebra. This needs to be studied on a case by case basis. Once this is 
achieved, the form factors for that operator can be determined. Concluding remarks are 
given in section 8. 

2 Model Hamiltonians 

Here we present, and give a description of, three models for Bose-Einstein condensates 
and the reduced BCS model. Our main objective is to establish that each model is 
exactly solvable through the algebraic Bethe ansatz. Throughout, there are no constraints 
imposed on the coupling parameters for all models other than they are real, which is to 
ensure hermiticity. 
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2.1 A model for two Josephson coupled Bose-Einstein conden- 
sates 

Consider the following general Hamiltonian describing Josephson tunneling between two 
coupled Bose-Einstein condensates 

H = U u N^ + U 12 N 1 N 2 + U^Ni+n 1 N 1 +n 2 N 2 

— y( a l a 2 + °2°l) C 1 ) 

where the operators Oj, aj, N{ = a\ai are associated with two Heisenberg algebras with 
relations 

[(hi a]} = [o-i, aj] = [a\, a ]} = 0- 

The Hilbert space of states is given by the infinite-dimensional Fock space spanned by 
the vectors 

\m,n) = (a t 1 ) m (4)" ' |0) , m, n = 0, 1, 2, oo. (2) 

The model describes Josephson tunneling between two condensates with tunneling strength 
Sj/2, the parameters U# are the amplitudes for S-wave scattering and /ij are chemical 
potentials. The Hamiltonian commutes with the total particle number N = N± + N 2 . 

The above Hamiltonian under the constraint U = Uu = U 22 = —U\ 2 /2 has been 
studied widely using techniques other than the exact solution |JSJ EH E3 ESI EH EE]- 
For this case it is useful to divide the parameter space into three regimes; viz. Rabi 
(U/Sj « N' 1 ), Josephson (A^ _1 « U/Sj « N) and Fock (N « U/Sj). In the Rabi 
and Josephson regions one expects coherent superposition of the two condensates to be 
possible whereas in the Fock region the two condensates will be, in some sense, localised. 
There is a correspondence between (P) and the motion of a pendulum [3H]- I n the Rabi 
and Josephson regions this motion is semiclassical, (i.e., the energy level spacings are of 
order less than N,) in contrast to the Fock case. For both the Fock and Josephson regimes 
the analogy corresponds to a pendulum with fixed length, while in the Rabi regime the 
length varies. An important problem is to study the behaviour in the crossover regimes, 
particularly between the Josephson and Fock regimes which are the most likely to occur in 
an experimental context [3B] ■ A reliable method to do this is through the exact solution. 
The motivation to extend the solution to the case where the couplings Uu, U 22 , U\ 2 for the 
S'-wave scattering terms can be chosen arbitrarily is for the description of a pair of Cooper 
pair boxes with capacitive coupling [IS!- I n the limit U 22 0, then (N 2 ) » (Ni), in 
which case the model can be considered as a single Cooper pair box coupled to a reservoir. 

2.2 A model for homo- atomic- molecular Bose-Einstein conden- 
sates 

Next we turn our attention to a two-mode model for an atomic-molecular Bose-Einstein 
condensate with identical atoms. The Hamiltonian takes the form 

H = U aa N 2 a +U ac N a N c + U cc N 2 c + fi a N a + /i c iV c 

+Q(a) a) c + aa) (3) 

which acts on a basis of Fock states analogous to (J2J. Here, a) is the creation operator for 
an atomic mode while creates a molecular mode. The parameters Uj again describe 
S'-wave scattering, /ij are chemical potentials and Q is the amplitude for interconversion 
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of atoms and molecules. The Hamiltonian commutes with the total atom number N = 
N a + 2N C . 

In the limit U aa = U ac = U cc = this model was studied in j^Ij, and analysed 
numerically in [SS| based on the Bethe ansatz solution. However, in order to compare 
with experimental results, in which the S'-wave scatterings are significant, one needs to 
analyse (j2J) in its full generality. Estimates for the S'-wave scattering parameters in the 
case of 87 Rb are given in [S3] . 

2.3 A model for hetero- atomic- molecular Bose-Einstein conden- 
sates 

The previous model can be extended to describe an atomic-molecular Bose-Einstein con- 
densate with two distinct species of atoms, denoted a and b, which can combine to produce 
a molecule c. For this case the Hamiltonian takes the form 

H = U aa N 2 a +U bb N* + U cc N 2 c +U ab N a N b + U ac N a N c + U bc N b N c 

+fi a N a + /j b N b + /i c iV c + ^(aWc + c ] ba) (4) 

which commutes with 1 = N a — N b and the total atom number N = N a + N b + 2N C . Here 
the model acts on the Fock space spanned by the vectors 

\l,m,n) = (o t ) , (6 t ) m (c t ) n |0) . 

Let us point out that in the limit U aa = U bb = U cc = U ab = U ac = U bc = 0, equation (J3J) 
is the Hamiltonian studied in (HH IHH modelling second harmonic generation in quantum 
optics. Non-zero values of these parameters correspond to a Kerr effect. 

2.4 The reduced BCS model 

The physical properties of a metallic nanograin with pairing interactions are described by 
the reduced BCS Hamiltonian jlU] 

c c 

H = J2 e 3 n i - 9 Yl 4+4- c i- c i+- ( 5 ) 

j=l j,k=l 

Above, j = 1, ...,£ labels a shell of doubly degenerate single particle energy levels with 
energies €j and rij = c^ + Cj + + cj_Cj_ is the fermion number operator for level j. The 
operators Cj±, cj ± are the annihilation and creation operators for the fermions at level j. 
The labels ± refer to time reversed states. 

One of the features of the Hamiltonian (J3J) is the blocking effect. For any unpaired 
fermion at level j the action of the pairing interaction is zero since only paired fermions are 
scattered. This means that the Hilbert space can be decoupled into a product of paired 
and unpaired fermion states in which the action of the Hamiltonian on the space for the 
unpaired fermions is automatically diagonal in the natural basis. In view of the blocking 
effect, it is convenient to introduce hard-core boson operators bj = Cj_Cj + , 6] = c] + c]_ 
which satisfy the relations 

(b]Y = 0, fo, bl) = Ml - 2b]b 3 ) [bj, b k ] = lb], b\] = (6) 

on the space excluding single particle states. In this setting the hard-core boson operators 
realise the su{2) algebra in the pseudo-spin representation, which will be utilised below. 
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The original approach of Bardeen, Cooper and Schrieffer to describe the phe- 
nomenon of superconductivity in a bulk system was to employ a mean-field theory using 
a variational wavefunction for the ground state 

c 

|*) = JJW + ^I)|0> (7) 
t=i 

which has an undetermined number of electrons. The expectation value for the number 
operator is then fixed by means of a chemical potential term /x; i.e. the grand canonical 
ensemble is used. One of the predictions of the BCS theory is that the number of Cooper 
pairs in the ground state of the system is given by the ratio A/d where A is the BCS 
"bulk gap" and d is the mean level spacing for the single electron energies. For nanoscale 
systems, this ratio is of the order of unity, in seeming contradiction with the experimental 
results discussed above. The explanation for this is that the mean-field approach is 
inappropriate in this instance, as previously indicated. 

3 Quantum inverse scattering method 

The essential motivation for the quantum inverse scattering method is the construction of 
a family of commuting matrices, known as transfer matrices. That is, we wish to construct 
an operator t(u), where u G C is called the spectral parameter, acting on some vector 
space, which represents the Hilbert space of physical states. Further we seek that 

[t(u), t(v)} = \/u,v E C. (8) 

There are two significant consequences of (|Sji. The first is that t(u) may be diagonalised 
independently of u, that is the eigenvectors of t(u) do not depend on u. This is the feature 
which makes the Bethe ansatz approach viable. Secondly, t(u) commutes with all of its 
derivatives, or more formally, taking the series expansion 

oo 

t(u) = ^ uk 

k=—oo 

it follows that 

Thus for any Hamiltonian which is expressible as a function of the operators only, each 
£fc corresponds to an operator representing a constant of the motion, since it will commute 
with the Hamiltonian. When the number of independent conserved quantities is equal to 
the number of degrees of freedom of the system, the model is said to be integrable. 

Let V denote some fixed vector space of finite-dimension n. The theory of exactly 
solvable quantum systems in this setting begins with an invertible operator, depending 
on the spectral parameter u, 

R(u) E End(V ® V) 

called the i?-matrix. Here "End" refers to the space of endomorphisms (square matrices), 
so R(u) is effectively an n 2 x n 2 matrix whose entries are scalar functions of u. From 
the i?-matrix we define the Yang-Baxter algebra, denoted Y, which is generated by the 
monodromy matrix T(u), whose entries are elements of Y 

R 12 (u - «)Ti(u)T 2 (i;) = T^T^R^u - v). (9) 
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The above equation acts in the three-fold space End(V r <S> V) <E> Y and the subscripts refer 
to the components of End(V ® V). In terms of the elementary matrices e*, which have 1 
in the position and zeros elsewhere, we may write 

n 

R(u)= i#(«)ej®e?, 

i,j,k,l=l 
n 

r(i») = ^e}®2j(«). 

Then 

^(w) = Yl R t^ 4 ® e ^ ® /; 

i ,j,k,l=l 
n 

n 

T 2 {u) = ^ J®e}®2j(u) 
where / is the identity operator. In component form we may write 

n n 

Y RfM ~ v)Ti(u)T l q (v) = T 3 k (v)Ti(n)R%(u - v) (10) 

3,1=1 j,l=l 

so the R l ^{u) give the structure constants of the algebra. Note that Y is actually an 
infinite-dimensional algebra, a basis for which {Tj[fc]} is obtained by taking the series 
expansions 

oo 
k=—oo 

Imposing that Y is an associative algebra leads, through repeated use of Q, to the 
following equation in End (V <E> V (g> V) (g> Y 

= (T 1 (u)T 2 (v))T 3 (w) 

= R£{u - v) (T 2 (v)Ti(u)) T 3 (w)R 12 (u - v) 

= = Ri 2 {u- v)R^(u - w)R^(v - w^w^v^u) 

x R 23 (v -w)R 13 {u-w)R 12 {u-v). (11) 

Here Rj^{u) denotes the matrix in End(V" ® V ® V) acting non-trivially on the j-th and 
k-th spaces and as the identity on the remaining space. In a similar way one may deduce 
that 

= T l {u) (T 2 (v)T 3 (w)) 

= R£(v - w)R u \u - w)R u \u - u)T 3 («;)T 2 (t;)Ti(u) 

x R X2 {u - v)R 13 (u - w)R 23 {v - w). (12) 
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A sufficient condition for and (|12|) to be equivalent is that the i?-matrix satisfies the 
Yang-Baxter equation acting in End(V <8> V ® V) 

Ru(u - v)R n (u - w)R 2 3(v — w) = R 2 z{v - w)R 13 (u - w)R X2 {u - v). (13) 

The above shows that in this algebraic setting the Yang-Baxter equation arises as a natural 
way to impose associativity of the Yang-Baxter algebra Y . It also appears in many other 
contexts, such as classical two-dimensional statistical mechanics [7] , knot theory [HJ [73] 
and scattering theory [73] . 



12 3 12 3 




3 2 1 3 2 1 

Figure 1: Graphical representation of the Yang-Baxter equation. In the context of 
knot theory, Rij(pd) denotes the crossing of string i over string j. The Yang-Baxter 
equation states that the above two combinations of crossings are topologically equivalent. 
In scattering theory the matrix elements of R%j{u — v) give the amplitudes for the two- 
body scattering of particles labelled i and j, with rapidity variables u and v respectively. 
The fact that the scattering depends only on the difference u — v is a consequence of 
Lorentz invariance. The Yang-Baxter equation is a statement of equivalence for the two 
factorisations of three-body scattering in terms of two-body scattering. In classical two- 
dimensional statistical mechanics the matrix elements of Rij(u) give the allowed vertex 
weights at the lattice site labelled by (i, j). In this instance u can be parameterised 
in terms of the energy levels and temperature. The Yang-Baxter equation ensures the 
commutativity of the row-to-row transfer matrix, from which the partition function is 
constructed. 

Here, we will only concern ourselves with the gl(2) invariant i?-matrix, which has the 
form E] 

R(u) = —^—(u.I^I + riP) 

U + 7] 
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/ 1 \ 

b(u) c{u) 
c(u) b{u) 

\ oiy 

with b(u) = u/(u + r]), c(u) = r]/(u + f]) and 77 is an arbitrary complex parameter 
P is the permutation operator which satisfies 

P(x ®y)=y®x Vx, y 6 V. 

The i?-matrix is gl(2) invariant in that 

[R(u), 0®fl] = 

where g is any 2x2 matrix. 

For this case the Yang-Baxter algebra, denoted Y[gl{2)} has four elements 



T{u) 



Tl{u) T%{u) 
T 2 Au) Ti(u) 



For clarity and convenience we adopt the notation 

A{u)=Tl{u) ) B(u)=T}(u), C(u)=T?(u) 
The full set of algebraic relations governed by (JHJ) are 



D{u)=Tl{u). 



[A(u), A{v)\ 
[B(u), B(v)\ 

[A(u), D(v)) 



A(u)B(v 
A{u)C(v 
D(u)B(v 
D(u)C(v 
B(u)A(v 
B{u)D{v 
C{u)A{v 
C{u)D(v 
[B(u), C(v)} 



[D(u), D(v)] = 

[C(u), C{v)\ = 0, 
V 



u — V 

7] 



(C(v)B(u) - C(u)B(v)) 



u — V 

U — V — 7] 



(B(u)C(v) - B(v)C(u)) 
V 





u 


— V 


u 




V + 7] 




u 


— V 


u 




V + 7] 




u 


— V 


u 




V — 7] 




u 


— V 


u 




V — 7] 




u 


— V 


u 




v + rj 
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— V 


a 




V + 7] 




II 


— V 


a 




V — 7] 



u — V 

7] 



B(v)A(u) + 
C{v)A{u) - 
B(v)D(u) - 
C(v)D(u) + 
A(v)B(u) + 
D(v)B(u) - 
A{v)C{u) - 
D(v)C(u) + 



u 




V 




V 




u 




— 

V 




V 




u 




V 




V 




u 




V 




V 




u 




V 




7] 




II 




V 




V 




a 




V 




V 




u 




V 



B(u)A(v), 
C(u)A(v), 
■B{u)D{v), 
■C(u)D(v), 
A(u)B(v), 
■D(u)B(v), 
A{u)C{v), 
■D(u)C(v), 



u — V 
V 

u — V 



(A(u)D(v) - A(v)D(u)) 
(D(v)A(u) - D(u)A(v)) . 



12 



Next suppose that we have a realisation of y[gZ(2)] acting on some vector space W, 
which we denote ix : Y[g/(2)] — > End W. It is usual to refer to V as the auxiliary space 
and W as the physical space. Note that as y[g/(2)] is a quadratic algebra, any realisation 
can be multiplied by an overall scaling factor and still satisfy the relations (fT7|). For later 
convenience we set 

L{u) = vr (T(u)) G End (V ® W) 
which we refer to as an L-operator. Defining the transfer matrix through 

t{u) = ix (tr (T(u))) = tx (A(u) + D{u)) G End W (18) 

it follows from © that the transfer matrices commute for different values of the spectral 
parameter; viz. equation (|SJ) is satisfied. 

An important property of the Yang-Baxter algebra is that it has a co-multiplication 
structure which allows us to build tensor product realisations. In particular, given two L- 
operators L u G End (V®U) and L w G End (V®W), then L = L U L W G End (V&U&W) 
is also an L-operator as can be seen from 

Rn(u - v)L 1 (u)L 2 {v) = R 12 (u-v)L^(u)LY(u)L^(v)Lf(v) 

= R 12 (u-v)L^(u)L^(v)LY(u)LY(v) 

= L^(v)L^u)R 12 (u~v)LY(u)Lf(v) 

= LV(v)LV(u)LY(v)LY(u)R 12 (u-v) 

= L^(v)L 2 v (v)L^(u)LY(u)R 12 (u-v) 

= L 2 (v)Li(u)Ri 2 (u-v). 

Furthermore, if L(u) is an L-operator then so is L(u + a) for any a, since the i?-matrix 
depends only on the difference of the spectral parameters. This property will prove 
important in all constructions below. 

3.1 The quasi-classical limit 

The i?-matrix (|14|) has the property 

lim R(u) = I ®I 

which is known as the quasi-classical property. For any such i?-matrix it is appropriate 
to write 

R(u) = I®I + nK(v) + o(r] 2 ) 
T ;(u) = Sil + rjTji^ + oirj 2 ) 

and substitute into (fTUj). Equating the second order terms in r\ yields the following 
relations 

n 

+K%(u - v)T*{v) - K%(u - v)7f{v)) 

which we take to be the defining relations for the algebra denoted y, to be called the 
Gaudin algebra. Gaudin used the quasi-classical limit to define classes of integrable spin 
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chain Hamiltonians with long-range interactions [7J3 EH] • The algebraic approach which 
we follow here is due to Sklyanin pTT) 177] . Observe that in the quasi-classical limit y is an 
infinite-dimensional Lie algebra in contrast to the quadratic algebra structure of Y. Any 
realisation of Y which admits the quasi-classical limit provides a realisation of y. 
For the case of the gl{2) invariant i?-matrix (|14|) let us write 

A(u) = I + r)A(u) + o(?7 2 ), 

B(u) = 7]B(u)+o(7] 2 ), 

C{u) = r]C(u)+o(r] 2 ), 

D(u) = I + r]V(u) + o(t] 2 ). 

From (|17|) we determine that the full relations for the algebra y[gl(2)} are 



[A{u) : 


■ A(v)] 


= [B{u) 


, B(v)} = 








[C(u) 


, C(v)} 


= P{u) 


, Viv)\ = 








[A(u)., 


V{v)\ 


= o, 










[B(u) 


, C(v)] 


A(u) 


- A{v) 4 


-V(v)-V(u) 








u - 


• 

- V 






[A(u) 


, B(v)] 


B{u) 
u 


-B(y) 

J 

— V 


[A(u), C(v)} = 


C(v) 
u 


-C{u) 
— V 


[D(u) 


, B{v)\ 


B(v) 
u 


-B{u) 

— V 


[V(u), C(v)} = 


C(u) 
u 


-C(v) 
— V 



3.2 Examples of realisation of the Yang-Baxter algebra 

In order to construct a specific model, we must address the question of determining a 
realisation of the Yang-Baxter algebra. Here we present several examples. The first 
realisation comes from the i?-matrix itself, since it is apparent by setting w = in (|TH|) 
that we can make the identification L(u) = R(u) such that a realisation of (jHJ) is obtained. 
This realisation satisfies the quasi-classical property, and is that used in the construction 
of the Heisenberg model [TU1 [TTJ H21 [T51 E]. A second realisation is given by L(u) = q 
(c- number realisation), where q is an arbitrary 2x2 matrix whose entries do not depend 
on u (although can depend on rj). This follows from the fact that (|T3|) holds for any such 
0- 

There is a realisation in terms of canonical boson operators a, with the relation 
[a, of] = 1 which reads [75] 

i»=( (1+ "; , a / * + " 2iV t ) (w) 

where N = a^a. There also exists a realisation in terms of the su(2) Lie algebra with 
generators S z and S ± [101 ELU EE2 113 Ej , 

-3/ \ _ 1 ( ul + r]S z r]S~ 



subject to the commutation relations 

[^,5 ± ]=±5 ± , [S + ,S-] = 2S Z . (21) 
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When the su(2) algebra takes the spin 1/2 representation the resulting L-operator is 
equivalent to that given by the i?-matrix. Another is realised in terms of the su(l,l) 
generators K z and K± [731 ISO] . 

- ; ( '"^f JT^. ) ■ m 

with the commutation relations 

[K z , K ± ] = ±K ± , [K + , K~\ = -2K Z . (23) 

Each of the realisations L a (u), L s {u) and L K {u) satisfy the quasi-classical property, and 
thus affords a realisation of the Gaudin algebra. 

The discerning reader may notice that L s (u) and L K (u) are in fact equivalent, which 
results from the Lie algebra homomorphism T : su(2) — > su(l, 1) defined by 

T(S Z ) = K z , T(S+) = -K + , T(S~) = K~ 

such that (|21|) is mapped to (|23p. For convenience we make the distinction between these 
two L-operators as the transformation T is non-unitary. This permits us to avoid the use 
of non-unitary realisations of the su(2) algebra below. (Although, as will be seen, the 
realisations of the Yang-Baxter algebra may not be unitary.) 



3.3 Z-graded structure and Z-graded realisations of the Yang- 
Baxter algebra 

The Yang-Baxter algebra Y[g/(2)] carries a Z-graded structure that can be exploited in 
the construction of models of Bose-Einstein condensates, as explained in [Hnj, which we 
now recount. We introduce an auxiliary operator Z, called the grading operator, satisfying 
the relations 

[Z,X(u)]=p{X(u)}.X(u), (24) 

where X = A, B, C or D and 

p{A(u)} = p{D(u)} = 0, p{B(u)} = 1, p{C(u)} = -1. 

We call p{X{u)} G Z the gradation of X(u), and extend the gradation operation to the 
entire algebra by the requirement 

P {6.(f)} = p{e}+p{(f)} ve,<f>eA. 

This definition for the grading operator is consistent with the defining relations (|17|). 

Let us now define a class of realisations ofY[gl(2)] which we call Z-graded realisations. 
We say that a vector space W, equipped with an endomorphism z, is a Z-graded vector 
space, denoted (W,z), if it admits a decomposition into subspaces 

oo 

w = w j 

j=-oo 

such that 

zW j = J .W j , jez. 
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Note that some of the Wj may be trivial subspaces. Formally, the grading operator can 
be used to define the following projection operators 

fc — — OO ' 

such that 

P k Pj = 5 kj Pj, P k Wj = 5 kj W k . 
We say that a Z-graded vector space 



is equivalent to W if for some k G Z there exist a vector space isomorphism between 
Wj and Wj+fc for all j. This terminology is motivated by the fact that for a given (W, z) 
one can always generate another Z-graded space (W, z') through the mappings Wj — > 
Wj +k , z' — > z — kl for any k G Z. 

For a given Z-graded W we say that 7r : y[pZ(2)] — > End IV provides a Z-graded 
realisation of F [g/(2)] if ir(Z) = z and the relations (fTTll24j) are preserved. In such a case 
we can write 



n(X(u)) = X ^ 



3 



and the matrices X(u,j) satisfy 

X(u,j)W k = foij^k. 

More specifically, this means that for \ipj) G Wj we have 

ir{X{u)Y{v)) \^j) = X(u,j+p{Y(u)})Y(v,j) |^> • 

In view of the equivalence of Z-graded vector spaces defined above, there can also exist 
equivalent realisations. We can define a realisation tt' equivalent to n by specifying some 
k G Z such that 

tt'(Z) =ti(Z -kl) 

and for 



n'(X(u))= J2 X '( u ^) 



the matrices X'(u,j) are defined by 

X\u,j)=X(u,j + k) Wj G Z. 
For the Z-graded case we may express the transfer matrix as 

OO 

t( u )= t( u >j) 

j=-oo 

such that 

t{u,j)W k = ioij^k 
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and 

[t(u,j), t(v, k)] = Vj, k G Z, u, v e C. 

Since p{t(u)} = 0, the diagonalisation of t(u) is thus reduced to the diagonalisation of 
each of the matrices t(u,j) on the Z-graded component Wj, where we have 

[t(u,j),t(v,j)] = Vu,veC. 

We may restrict our attention to the case of t(u, 0), as each t(u,j) is equivalent to some 
t'(u, 0) through the use of equivalent realisations as introduced above. 

3.4 Examples of Z-graded realisations 

Next we give two non-trivial Z-graded realisations of the algebra Y[gl(2)]. One is ex- 
pressible in terms of two Heisenberg algebras with generators aj, a], i — 1, 2 and reads 
X(u,j)=X(u,j)P j with 



A{u,j)) 


= u 2 + r]uN + T) 2 NxN 2 - r?(iVi - 
-co 2 (N + jl) + 4a 1; 


N 2 )ou(N + jI) 


B(u,j) 


= (u + u(N + jI)+r]Ni)a 2 + r]' 


"V, 


C(u,j) 


= a\(u-u(N + jr)+r}N 2 )+r}' 




D(u,j) 






Z 


= k.I-N. 





Above, k is an arbitrary scalar, Pj are the projections defined by (J25)) . iVj = ajcij, iV = 
N\ + A^2 and u(x) is an arbitrary polynomial function of x. The operators act on the 
Fock space spanned by the basis vectors given by (J2J). Note that in the case when uj(x) 
is constant, the above realisation reduces to that discussed in jlHl EDI E2 E3 IMj and is 
factorizable into two local realisations of the Yang-Baxter algebra expressible in terms of 
the two Heisenberg algebras; viz. 

L{u) = rT 2 L\{u - r]- 1 + u)L a 2 {u - r]' 1 - u). 

It is important to note that for generic u(x) no such factorisation exists. 

Another Z-graded realisation of the Yang-Baxter algebra is X(u,j) = X(u,j)Pj with 

A(u,j) = -r]u 2 + u(l-r] 2 (K z + N c ) -rju(K z + N c + jl)) 

+r]K z - rfK z uj(K z + N C + jl) - r] 3 N c K z + r] 2 cK + , 
B(u,j) = T](l - T]U-T]U(K Z + N c + jl) - 7] 2 N C )K' - r]c(u- r]K z ), 
C(u,j) = r]c\u + r]K z ) -r]K + , 
D{u,j) = u-r]K z + r] 2 c l K~, 

Z = k.I - K z - N c + re. (27) 

Above, k is again arbitrary, the operators c, form a Heisenberg algebra, with iV c = c^c, 
and as before the operators K z , K + , K~ satisfy the su(l, 1) relations (|2^jl . It is assumed 
that the su(l, 1) operators are realised in terms of an irreducible representations of lowest 
weight re. As in the previous example, u(x) is an arbitrary polynomial function of x 
and the above realisation is factorisable only in the case when u(x) is constant. In this 
instance we have 

L{u) = uQL a (u - rT 1 + u)L K {u) 

where g = diag(— 1, 1). 
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4 Algebraic Bethe ansatz method of solution 

A key step in successfully applying the algebraic Bethe ansatz approach is finding a 
suitable pseudovacuum state, \x), which has the properties 

A(u)\0) = a(u)\ X ), 
B(u)\0) = 0, 
C(u)\0) ? 0, 
D(u)\0) = d(u)\ X ) 

where a(u) and d(u) are scalar functions. Note that for ease of notation throughout we 
will omit the symbol 7r denoting the realisation of Y[gl(2)]. Next choose the Bethe state 

M 

\^) = \v 1 ,...,v M ) = l[C(v i )\x). (28) 

i=l 

Note that because [C{u), C(v)] = 0, the ordering is not important in the product of (J25j) . 
The approach of the algebraic Bethe ansatz is to use the relations (|17j) to determine the 
action of t(u) on \v). First let us consider the action of A(u) on \v); viz. 



where 
Now 



A{u) \v) = A{u)C{ui) \vi) 
Wi) = \vt, ■ ■, Vi-x, v i+ x, v M ) . 



A(u)\v) = U Vt + V C(v i )A(u) \v t ) - -^—C(u)A( Vi ) \v t ) 

U — Vi U — Vi 

" ' ^^-CivAAiuMvj) \v tJ ) - -Jt—C(u)A( Vi )C( Vj ) \v t ,) 



U — Vi U — Vi 

C( Vi )C( Vj )A(u) |%) 



U — Vi + 7]\ ( U — Vj + T] 



U — Vi I \ U — Vj 



U — Vi + Tj 



U — Vi / \ U — Vj 
V \ f Vi - Vj + 7] 



U — Vi I \ Vi~ Vj 



C{vi)C{u)A{vj) \vij) 
C(u)C(vj)A(vi) \vij) 



U — Vi l \ Vi — Vj 



+ [ - 1 - ) ( — — ) C(u)C( Vi )A( Vj ) \v i5 ) , (29) 



where, 



1%) = \ V U - - - > Vi-i,V i+ x, Vj-x, v j+1 , •■■ , v M ) ■ 

Proceeding further we find the general form 

A{um = (n^^)(n^))^)ix) 

— ^- (f[ V -^l±A C(u) (f[C(v k )) A{vi) \ X ) 
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+ other linearly independent terms 



M 



a\u) 



it fj 

v«=l 



7 7 a (^i) ( T\ v i - v j + V 



J. J. 71. — 71. / 



AMP) = .wn^H 



W - 1>i \ f » - U, 

+ other linearly independent terms. 

Above, each of the other linearly independent terms is a vector of the form C(u) \vj) , j 
multiplied by some scalar. There are no other possibilities. To determine what the co- 
efficients are we note that the above equation is valid for any choice of i. Hence we 
conclude that 

M . s 

U -Vi + T) 

U-Vi 

f—' U - Vi \ 1 Vi- Vj I 

1=1 ■> / 

The co-efficients of the terms C(u) \vi) are called unwanted terms for reasons that will 
soon become apparent. 

We now perform the same procedure for D(u); 

D{u)\v) = D{u)C{v % )\vi) 

U ~ Vi ~ T] ) C{v,^D{u) \vi) + — — C{u)D{vi) \v t ) 

U-Vi J U-Vi 
M 



v i=l 



» x x U — Vi I 



+E 2*!0fn«izaz3| C( „ )B >. (30) 



. IV: Vi- vj 

1 = 1 \J7=* 



The final result for the action of the transfer matrix is 

t(u)\v) = (A(u) + D(u))\v) 
= A(u, v) \v) 



where 



+ f;!*ifn^i^l N ) C ( M )| i r i > (3i) 

*-fu-Vi \*J: Vi- Vj I 

i=l J / 

M M 

A(u, 0) = a(«) TT U ' Vl + V + d{u) TT ^"^'^ (32) 

- LJ - U — Vi - LJ - u — Vi 

1=1 i=l 

The above shows that becomes an eigenstate of the transfer matrix with eigenvalue 
(jS2)) whenever the unwanted terms cancel. This occurs when the Bethe ansatz equations 

t \ M 

Ci{Vi) t — r Vi — Vj — 77 

-tV^- = II — "-T 1 , i=l,...,M. 33 

a(%) j~ Vi — Vj + r] 
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are satisfied. Throughout we adopt the notation 



{v^ = {v x ,v 2 , ■ ■ ■,%} 



for such a solution. 

Note that in the derivation of the Bethe ansatz equations it is required that Vi 7^ 
Vj j. This is a result of the Pauli Principle for Bethe ansatz solvable models, as de- 
veloped in [HI] for the one-dimensional Bose gas with delta-function interactions. For 
generic functions a(u), d(u), essentially the same argument as [EI] can be applied to draw 
the same conclusion. To give an indication why this is the case, consider (129)) in the limit 



This equation shows that new types of unwanted terms occur which depend on the deriva- 
tives of the elements of Y[gl(2)}, and this leads to an overdetermined system of equations 
which do not admit a solution. Another viewpoint is to note that (up to an overall scaling 
factor, and under suitable assumptions for the forms of a(u) and d(u)) the eigenvalues 
A(u, v) are analytic functions of u. Assuming that the poles in ()32|) are simple then the 
Bethe ansatz equations (j3*3|) are equivalent to the statement that the residue vanishes at 
each pole; i.e., 



leads to For non-simple poles however there are additional Bethe ansatz equations, 
which cannot be satisfied. To illustrate this, we consider the simplest case where, for all 
values of 77, all poles are simple except for one, say at u = Vj, which is second order. In 
such a case, let us write 




lim {u — Vi)A(u, v) = 



(u - Vj + rf) 

(U - Vj) 2 



2 M 



(ll-V k + Tj) 
(u - V k ) 




A(w, v) = a{u) 



n 



+ d(u) 



Analyticity of A(u, v) requires 




leading to the Bethe ansatz equations 




M 




Vj -Vk-y 

Vj -V k + 7] 



In particular note that 




-1. 



(34) 
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Furthermore 



= lim (u — Vj)A(u, v) 

U->Vj 

= lim — [(it — Vj) 2 A(u, vj\ 

/ M M 

= lim a'{u){u - Vj + r]) 2 JJ(u - v k +rj) + 2a{u) \\{u -v k + r)) 

u ~* Vj \ k^j k=i 

M M 

+ a(u){u - Vj + v) 2 ^2 II ( n ~~ Vk + 

M M 

+ d'(u)(u - Vj - r]) 2 Y\(u -v k -T]) + 2d(u) \\{u - v k - rj) 

k=ij k=l 
MM \ 

+d(u)(u - vj - v) 2 ^2 Y[(u - v k - 7]) \ 

¥3 M« / 

M M 

= 2r]a(vj) Y[(vj -v k + r))- 2r]d(vj) Y[( v 3 ~ v k ~ v) + (v 2 )- 

k^j k^j 

The above implies that 

lim -p( = 1 

v^o dyVj) 

in contradiction with (|34|) . proving the claim that one cannot have a single second order 
pole. The extension to more complicated non-simple pole structures, though tedious, is 
straightforward. 

Finally, the following identity is useful and is easily derived from the Bethe ansatz 
equations 

M M 

n a K>=n^)- ( 3s ) 

i=l i=l 

4.1 Extension to Z-graded realisations 

In order to formulate the algebraic Bethe ansatz solution for the class of Z-graded realisa- 
tions, we begin with the observation from (jl7j) that the following relations hold (amongst 
others) : 

[A(u,j),A(v,j)] = [D(u,j),D(v,j)] = 0, 
B{u,j)B{v,j-l) = B{v,j)B{u,j-l), 
C(u,j)C(v,j + l) = C(v,j)C(u,j + l), 

A(u,j)C(v,j + l) = U ~ V + V C(v,j + l)A(u,j + l) 

u — V 

-C(u,j + l)A(vJ + l), 



u — V 



D(u,j)C(v,j + l) = U V V C(v,j + l)D(u,j + l) 

u — V 



.^!—C(u,j + l)D(v,j + l), (36) 
u — V 
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Again, we assume the existence of a pseudovacuum vector \x) £ Wk such that 

A{u,k)\x) = a(u,k)\x) 

B{u,k)\ X ) = 

C(u,k)\ X ) ± 

D(u,k)\ x ) = d(u,k)\ X ). 

In particular, for the realisation (|2*U|) the Fock vacuum serves as the pseudovacuum. In 
the case of (JTTj) we choose the pseudovacuum to be the tensor product of the Fock vacuum 
with the su(l, 1) lowest weight state of weight n. 

The above implies that \x) is a maximal weight vector with respect to Z. Without 
loss of generality we can choose k = M, due to the equivalence of realisations discussed 
earlier, and look for Bethe states defined by 

\v) = C(v 1: l)C(v 2 ,2)---C(v M ,M)\x). (37) 

It is easy to check that the above Bethe state is symmetric with respect to interchange of 
the variables Vi, a feature which plays a crucial role below. In particular, this means that 
we may write 

\v) = C(v h 1) \vi) 

= CK 2)|"%> 
= C{v v l)C{v u 2)\v tJ ) 

where now 

= C(v h 2)---C( Vi „ h i)C(v i+ i,i + l)---C(v M ,M)\x) } 
|%) = C(v 1 ,3)---C(v i ^ 1 ,i + l)C(v i+1 ,i + 2) 

x • • • C(^_i, j)C(v j+1 , j + 1) • • • C(v M , M) \ x ) . 

Acting A(u, 0) and D(u, 0) on the Bethe state we have, by following the general pro- 
cedure detailed above, 

A(u,0)\v) = a{u,M) \Y\ U ~ Vt + V )W} 

rja(vi,M) lYTVi-vj +r] 



_y <]^ll^± TT C(«,l)|«i), 
D(u, 0) \v) = d(u,M) \ f[ U ~ V ^~ r] ) \v) 



vi=l 



T)d(Vi, M) ( t— r Vi Vj T] 



+ y n C (u, i) %) . 

^— ' II 11: \ J-J- 11: 11: I 



U — Vj \ Vj — Vj 

i=l 1 



Requiring \v) to be an eigenstate of t(u, 0) leads to the Bethe ansatz equations 



^j = n^^, i = l,-,M (38) 
d(Vi, M) j~ Vi -Vj+7] 



and the corresponding eigenvalue of the matrix t(u, 0) is 



M M 



A(u, 0, v) = a(u, M) TT U Vi + V + d(u, M) TT U Vi V . (39) 

J. J. ni — 71. J. J. ni — 71. 



U Vi ^-^ U-Vi 

1=1 1=1 
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5 Scalar products of states 

Recall that in the usual algebraic Bethe ansatz for the algebra Y[g/(2)] there is a formula 
originally due to Slavnov [32] (see also O I2H]) for the wavefunction scalar products, 
which is 



S(w : v) = {w\v) 
= {v\w) 



AetF (40) 



Hk>l( V k - V^IliKjiWi - Wj) 

with the entries of the M x M matrix F given by 

d(w) ( M M \ 

F v = 7 — Z~\ a ( v i} ~ w k + v) - d{vj) Y[{vj - w k - r))\ (41) 

^ Wi > \ k+i k^i J 

and {u\ is the left vector defined by 

{u\ = ( X \B{u M )---B{ Ul ) 

for any choice of Above, {wi} provide a solution to the Bethe ansatz equations (JHHj) 
and the parameters {vj} are arbitrary. In using the Slavnov formula it is assumed that 
the pseudovacuum state has norm equal to one. Defining 

G = FT 

where T is a diagonal matrix with entries 



Uk=l( V j ~ W k) 

the Slavnov formula may be expressed in the equivalent form 



S{w:v) = — M - -detG (42) 

Uk>i( v k ~ vi) [[^{w, - Wj ) 



with 



G 



rjd{wi){vj — w 



13 (Vj - Wi) 2 

<*»n^- %) 5^)' (43) 

We will find it convenient to use both forms ()40|42|) of the Slavnov formula. 

The Yang-Baxter algebra F[gZ(2)] admits a conjugation operation f : Y[g/(2)] — > 
Y[gl(2)} defined by 

A(«) f = A(u*), B( U y = C(u*), C( U y = B(u*), D( u y = D(u*) 

and extended to all of Y[gl{2)] through 

(0.</>)t = 0tjt ; ye, 4>eY[ g i(2)} 
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such that the defining relations (|17j) are preserved. Above, * is used to denote complex 
conjugation. Consequently the right vector 

\v} = {v\ 

= B(viy---B(v* M y\ x } 

is also an eigenvector of the transfer matrix whenever the Bethe ansatz equations for 
the parameters {v i} are satisfied. However, it is apparent that the Z-graded realisations 
I27|) we have introduced are not unitary, and generally 



(v\ = \v) 7^ {v\ . 

On the other hand, numerical analysis for the models ()1I3I4I5|) indicates that for fixed 
particle numbers, and generic values of the coupling parameters, the energy spectrum 
is free of degeneracies. This is presumably due to the fact that the only Lie algebra 
symmetries for these models are u(l) invariances corresponding to conservation of particle 
numbers, and the non-degenerate spectra are examples of Hund's non-crossing rule [E21 
We also assert that for a given {fj}, satisfying the Bethe ansatz equations (|5^j). this 
set of parameters is equivalent to {v *}; i.e., v * = Vj for some j = 1, • • • , M. It is clear 
that for {vi} satisfying so does {v*} (in all our examples, as the Hamiltonians are 
real, the functions a(u) and d(u) are real, as will be seen below). Since the eigenvalues of 
the Hamiltonian are real, we have 

E(v) = E*(v) = E(v*). 

Under the belief that the spectrum is multiplicity free, we then deduce {v^} ~ {v*}- 
Whenever this is the case, we can conclude that the eigenvectors are real and 

(v\ = C(v) {v\ (44) 

for some non-zero real- valued scalar C(^)- Thus the Slavnov formula can be invoked for 
the computation of form factors. Throughout we will always assume (pHj) to be the case, 
which implies that 

{w\v) = for v 7^ w (45) 

whenever {v{\ and {wj} both satisfy the Bethe ansatz equations. This result (J45j) can be 
proved directly, independent of (JSJ), as shown in [HOj. Note when w = v we need to take 
a limit for the diagonal entries of F to compute the square of the norm. This yields 

(M M 
a '( v i) Yl( V i ~ V k + V) ~ d'( v i) Y[( V i - V k-V) 
k=l k=l 

(M / x M 

_^_«_JJ ( ) (46) 

where the prime denotes the derivative. 

From the Slavnov formula the matrix elements of the operators B(u), C(u) follow 
directly, as the set of parameters {v{\ are arbitrary. We will now derive an expression for 
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the form factors of the operator D(u), which will prove useful for later calculations. From 
the expression (JHUj) for the action of D(u) on an arbitrary Bethe vector we may deduce 
that for both {v^ and {wj} satisfying the Bethe ansatz equations 



{w\D(u)\v) 



M 



d( U ) n 



U — Vp — Tj 



M 



U-Vp 
p=l K 

M 



{w\v) 



Vg -Vj-T) 



U Vq 
q=l y 



v q - Vj 



{w\C{u)\vi) 



M 



d{u) n 



11 — Wp — Tj 



U-Wp 
p=l H 



S{w : v) 



M 



r]d{v q 



M 

n 



u-w k 



' x {u- v q )d(u) V £A u-w k -T] 



x S(w : Vi, ■ ■ -,v q _i,u,v, 



M 

n 



v q -Vj-ri 

Vg - Vj 



d{u) 



q— 1) ">■> u q+ll ' j 

M 



V M ) 



M 



Uk>i( v k - vi) n^,K - W. 



n 



U — Wp — Tj 



J J \ p= l U ~ W P 



lk>l\" K " l ) 1 H<j 

where the matrices F^ g > are defined by 

b pl - 

Fpq = Qpqi 

Q is the rank one matrix with elements 



detF + 0(w)^detF {9) 

9=1 



F p i for I ^ q 



(47) 



Qij 



r]d(wi)d(vj) 



(u - Wi){u - Wi - Tj) 



1 - 



and 



Q(u) 



a[u) 



M 

n 

vfc=l 



M 

n 



U-W k + T) 
U -W k -T] 



M 



Y[{vj -vi-rj) 



i=i 



u-w k 
u v k 



(46 



Using the fact that for v ^ w, det F — 0, whereas for v = w, Q(u) = 1, allows us to write 
{w\D{u)\v) 



®{u)d{u) 



M 

n 

M 

Uk>i( v k - vi) wTAwi - Wj) 



Q(u)d(u) 



U — Wi — 1] 



u — W,; 



U — Wi — Tj 



U — Wi 



M 



det F + det F 

\ 3=1 

det (F + Q). 



(?) 



(49) 



The last line above follows from the fact that if X is any M x M matrix and Y is any 
rank one M x M matrix then 



M 



det(X + Y) = det X + ^ det 
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where 



X$ = X kl for j ^ I. 

A similar result can be derived for the form factors of A(u). Alternatively, one can obtain 
them from (|32I49|) through 

{w\A(u)\v) = A(u,v) {w\v} - {w\D(u)\v} . (50) 

The Slavnov formula can be extended to include generic Z-graded realisations as noted 
in |66j . This is achieved by simply replacing a(u) and d(u) with a(u,M) and d(u,M) in 
(j41J43j) . All the results derived above also extend analogously to the Z-graded case. 

5.1 The quasi-classical limit 

Assuming that the quasi-classical limit exists and in particular 

a(u) = 1 + T)a(u) + o(r] 2 ) 

d(u) = l + rjd(u) + o(r] 2 ) 

it is straightforward to obtain the Slavnov formula in the quasi-classical limit. We obtain 
directly from ()40I42|) 

S(w:v) = (x\B(w M )---B(w 1 )C(v 1 )---C(v M )\x) 
= (x\B(v M )---B(v 1 )C(w 1 )---C(w M )\ X ) 

AV n — Wr,) 

■detg 



T-rM t-tM , \ 



Uk>l( V k - Vl)Ui<j(Wi ~ w j) 

— m 17 detjF (51) 



i<3 

where the entries of the M x M matrices Q and T are given by 



Qa = ( «(%■) - ' + E — ^— ) r Vj % ( 52 ) 

Above, the parameters {wj} are to satisfy the quasi-classical limit of the Bethe ansatz 
equations 

M 

a{ Wi )- t( Wi ) = y* , i = l,---,M (54) 

k^i 

while the set {vj} are arbitrary. 

Specialising to the case when {vi} = {wi} leads to the formula 

S(v : v) = det/C 
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where 



M 2 

Ka = a(Vi) - d'(Vi) - } 

K, i:j = — for i ^ j. 



We can perform a similar treatment to yield the quasi-classical limit of (|49|) . The 
terms in rf M give the scalar product of the states in the quasi-classical limit. The terms 
in r) +1 give not only the form factor for T>(u) but also the next order terms in the 
expansion of the scalar product, so some care needs to be taken in order to identify the 
appropriate terms. The result is 

{w\V(u)\v} - 



Uk>l( V h ~ Vl)Ui<j(Wi ~ Wj) 
M \ 

x[[d(u)-y \detJ r +y ' det F [q) I (55) 



M 1 \ M \ 

det^ + Vdet^) 



where J 7 ^ is defined in terms of T and Q in analogy with ()47|) and 

e*(«) = {u _ Wi? [*M -<(«)+ E ^ J ■ ( 56 ) 

6 Exact solution of the models 

6.1 Solution for the model of two Josephson coupled Bose-Einstein 
condensates 

It is an algebraic exercise to show that the Hamiltonian ((TJ is related with the matrix 
t(u, 0) = A(u, 0) + D(u, 0) obtained through dHJ) via 

# = — ^ p(0,0) -7]~ 2 + (aN + (3) 2 -7]aN-7]5N 2 ] 

where we have chosen u(N) = aN + f3 and the coupling constants are identified as 

2 2{U n + U 22 - U 12 ) 



V 



a 







a 



Noting that 



Hi + 1^2 

Un + U22 



N — 77 ^-(0, 0), 
au 



(57) 
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the above demonstrates that the Hamiltonian (|1J) is expressible solely in terms of the 
matrix t(u, 0) and its derivative. 

Since [H, N] = 0, the Hamiltonian is block diagonal in the Fock basis (0). Thus 
on a subspace of the Fock space with fixed particle number N, the diagonalisation of 
t(u, 0) is equivalent to the diagonalisation of t(u, 0) presented earlier in the Bethe ansatz 
framework. It is easily determined that for this case, the total particle number N = M 
and 

a(u,N) = u 2 - (aN + (3) 2 , 
d(u,N) = if 2 . 

From (jH8JH9j) we deduce the solution of (JTJ) for the energy spectrum to be 

N N 



E(v) = 
v ' 2 



j=l ' 1=1 



-tT 2 + (aiV + /3) 2 - W N - r]5N 2 ] (58) 
where the parameters {vi} are subject to the Bethe ansatz equations 

N 

V 2 (v 2 -(aN + ( 3) 2 )=H V *- V i- 7] . (59) 

' A j- V; — Vi + T] 

6.2 Solution for the model of homo-atomic-molecular Bose-Einstein 
condensates 

In terms of a realisation of the algebra su(l, 1) through 

K+ = W = £ K . = m 

one may establish the relation between the Hamiltonian (J3J) and the corresponding matrix 
t(u, 0) = A(it, 0) + D(u, 0) arising from the realisation (jTTji of the Yang-Baxter algebra is 

# = or + 5(iV/2 + 1/4) + 7 (iV/2 + 1/4) 2 + 2r]- 2 Qt(0,0), 

with 

-^(0, 0) = 2 - 77(77 + a)(JV/2 + 1/4) - v {3. 
au 

Above we have chosen 

u{K z + N c ) = a{K z + N c )+(3 
= a{N/2 + l/A)+ (3 
and the following identification has been made for the coupling constants 

4U aa + U cc - 2U ac 



V 
a 

P 
a 

5 

7 



2Vt 





2VL 


J 


2fj, c 




+ AU aa - U ac 






40 


u aa 


- 2/i a 
4 




2fi c 


-u ac 

2 




Ucc- 
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We deduce 



a(u, M) = (u + t)k) (1-rju-rj (a(M + k) + /?)) 
d(u, M) = u — r\K 

and by the same argument as in the previous example, we conclude that the exact solution 
for the energy spectrum of is determined by (|38I39|) which reduces to 



E(v) = ct + 6(M + k)+7(M + kY 



+2ti~ 1 kQ 



M M 



Vi-T] T-jVi+T] 



1=1 1=1 



(61) 



where the parameters {v satisfy the Bethe ansatz equations 



/ , \ M _ _ 
[1 - m - rj(a(M + (3)} ( ) = ]J Vi V > V . (62) 

3^ 

For the representation (|6T1|l of the su(l, 1) algebra there are two lowest weight vectors; 
viz. the Fock vacuum |0) and the one particle state |0). It follows from ()6())1 that the 
allowed values for k in ()61l62j) are k = 1/4, 3/4. This demonstrates that the solution of 
the model depends on whether the total particle number iV = 2M + 2k, — 1/2 is even or 
odd, the effects of which on the energy spectrum can be seen through numerical analysis 

(cf. |ESj). 



6.3 Solution for the model of hetero-atomic- molecular Bose- 
Einstein condensates 

In order to show the solvability of the model (J3J), we adopt the realisation of the su(l, 1) 
algebra given by 

K + = <Jb\ K- = ab, K z = iVa + ^ + 1 , (63) 

and observe that the operator X = N a — Nf, commutes with the su(l,l) algebra in this 
representation, hence taking a constant value in any irreducible representation. Due to the 
symmetry upon interchanging the labels a and b, we can assume without loss of generality 
that the eigenvalues of X are non-negative. In particular, note then that the lowest weight 
states for this realisation are of the form 

(a^) m 

|m) = ^=|0), m = 0,l,2,---,oo 



and K z \m) = (m/2 + 1/2) |m). We conclude that the lowest weight labels k can be taken 
from the set {1/2, 1, 3/2, ....} and the eigenvalue of X in the irreducible representation 
labelled by k is 2k — 1. 

For this case the relation between the Hamiltonian (0J) and the corresponding matrix 
t{u, 0) from (J27|) is 

H = a + 8(N/2 + 1/2) + \(N/2 + 1/2) 2 

+ pi + vi 1 + £X(iV/2 + 1/2) + r]- 2 nt(0, 0) (64) 
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with 



^(0, 0) = 2 - 77(77 + a)(N/2 + 1/2) - m I - V {3. 



Above we have chosen 

uj(K z + N c ) = a(K z + N c ) + 7 (2re - 1) + (3 
= a(N/2 + 1/2) +7J + /5 

and the coupling constants are related through 

U aa + U bb + U cc + ?7 afe — U ac — U bc 



V = 



a = 



7 



n 

Ucc — Uaa — Ubb — Uab 

2U aa + 2U bb + 2U ab - U ac - U bc + 2fi c - 2fi a - 2pt b 



2VL 



2U bb - 2U aa + U ac - U bc 



2VL 

U aa + U bb + U ab - 2jl a - 2fl b 



We find 



5 = 
A = 

P = 



2pc — Uac — Ubc 

2 ' 

U CC1 

U bb — U aa + fi a — Pb 



U aa + U bb — Uab 



U ac — U bc 



a(u, M) = (u + t)k) (l-rju-T] (a(M + re) + 7 (2re - 1) + (3)) , 
d(u, M) = (u- t}k) 

and the exact solution in this instance reads 

E(v) = a + 5(M + re) + A(M + re) 2 

+ p(2re - 1) + i/(2re - l) 2 + £(2re - 1)(M + re) 



+ 77" 1 refi 



M 



(1 - 77(a(M + re) + 7 (2re - 1) + f3)) \\ 



Vi-7] 



i=l 



Vi 



M 

n 

i=i 



where the parameters {t> j} satisfy the Bethe ansatz equations 



[1 -77^-77 (a(M + re) + 7 (2re - 1) + /?)] 



M 



Vi -Vj-T) 



n— 



Vi-T]KJ Vi -Vj+T) 



The total atom number is given by iV = 2M + 2re — 1. 
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6.4 Solution for the reduced BCS model 



As an alternative to the BCS mean-field approach, one can appeal to the exact solution 
of the Hamiltonian (jSJ) as given in [SHI dH|- Much later it was shown by Cambiaggio, 
Rivas and Saraceno [HI] that © is integrable in the sense that there exists a set of 
mutually commutative operators which commute with the Hamiltonian. Our aim here 
is to show that both of these features are consequences of the fact that the Hamiltonian 
(JSJ can be derived using the quantum inverse scattering method. This result, which 
was established in [6*T] l§2"] , will be proved below. Before doing so, let us remark that 
there have been several works on this problem, including generalisations [23 ESI EH] • In 
some cases these models can be obtained using trigonometric/hyperbolic versions of the 
Yang-Baxter algebra. While we will not go into details here, this generalisation from 
the procedure described below is straightforward and is simply a matter of using the 
trigonometric/hyperbolic analogue of (fTlj) from the outset. 

We use a c- number realisation g of the L-operator, defined by g = exp(— arja) with 
a = diag(l, —1), as well as (J2U|) to construct the transfer matrix 

t{u) = tr {q L s 0C {u - e c ) ■ - e t )) (66) 

which is an element of the £-fold tensor algebra of su(2). Here, tr denotes the trace 
taken over the auxiliary space, which for convenience is labelled by 0, while the tensor 
components of the physical space are labelled 1, • • - , C Defining 

Tj = lim tt^Uu) 

u-^ej rf 

for j = 1, 2, C, we may write in the quasi-classical limit Tj = Tj + o(i]) and it follows 
from the commutativity of the transfer matrices that [tj, Tfe] = 0, Vj, k. Explicitly, these 
operators read 

Ti = ~2aS* + J2 — (67) 
Mi 3 €k 

with 9 = S + <g> S~ + S~ ® S + + 2S Z ® S z . The set of operators (jBTj) . realised in terms 
of canonical fermion operators, are those obtained by Cambiaggio et al. 91j to establish 
the integrability of the reduced BCS model. They first appeared in the work of Sklyanin 
[77] in a general context, and are the Gaudin Hamiltonians [7o"[ I7B] in the presence of a 
non-uniform magnetic field. 

Next define a Hamiltonian through 



H 



c c c c 



a — - - 4a 3 ^ 3 2a 2 ^ 3 2a 

j=l j,k=l j=l j=l 

C 



5>^~E^ (69) 



a 

3=1 j,k=l 



where 



C = S + S- + S-S + + 2(S 



z\2 



is the Casimir invariant for the su{2) algebra. The Hamiltonian is universally integrable 
since it is clear that [H, Tj] = 0, Vj, irrespective of the realisations of the su{2) algebra 
in the tensor algebra. 
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In order to reproduce the Hamiltonian (j5J) we realise the su(2) generators through the 
(spin 1/2) hard-core boson representation ©; viz 

S-=b 3 , Sf = b], fi£ = i fo- - i) . (70) 

In this instance one obtains (jHJ) (up to the constant term — J2f=i e j) with g = 1/a. In- 
corporating higher spin representations of the su(2) algebra yields models which may be 
interpreted as coupled BCS systems EH EH1 Ell • Generally, we can define a represen- 
tation of su(2) through 

aeci o-en aeo. 

where the operators aj a , at may be either bosonic (+ sign in or fermionic (- sign 
in Sj). Above, a G Vt is a degeneracy label and a ^ Vt refers to the time-reversed state 
(i.e., the total degeneracy is twice the cardinality of Q). In this instance one recovers 
the pairing models discussed in (HSl EE1 EH] . (For the bosonic case it is convenient to 
replace L (u) with L K (u) since ()71|) is not unitary for bosons.) Because of these different 
possibilities we will derive the eigenvalues of the Hamiltonian (|69|) in a general setting. 

For each index k of the tensor algebra in which the transfer matrix acts, and ac- 
cordingly in (jnnj), suppose that we represent the su{2) algebra through the irreducible 
representation with lowest weight (or spin) — Note that we impose no restriction on 
the allowed values of s k in order to accommodate infinite-dimensional representations such 
as the bosonic case of (J7TJ). Choosing the pseudovacuum state to be the tensor product 
of lowest weight states gives 



a(u) 



c 

I \ TT u ~ e k - V s k 
exp(-ar/) ]_} 



k=i u ~ ek 

d(u) = exp(a?7) TT - — 6fc r l Sk 

u — efc 

fc=i 



and the eigenvalues of the transfer matrix ()66|) as 



c M 



u - e k - 7]s k -p-r u - Vj + r] 



J. J. Ol c, J. J. 



u - e k * * u Vj 

k=l j=l J 

C M 
/ \ TT U ~ € k + V S k -TT U ~ V j - 

+ exp(m?) I I I I J — 

7" \ u-e k . "t u Vj 

k=l 3=1 J 



The corresponding Bethe ansatz equations read 



exp 



C M 

v% — e k — r]s k T-r Vi — Vj 



, Vi - e k + 7]s k " Vi -Vj + rj 

k=l j=£i J 



The eigenvalues of the conserved operators (|57jl are obtained through the appropriate 
terms in the expansion of the transfer matrix eigenvalues in the parameter rj. This yields 
the following result for the eigenvalues Xj of Tj 



h= h«+E — -Er^r \; (72) 
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such that the parameters {vj} satisfy the quasi-classical limit of the Bethe ansatz equa- 
tions 

2a + V^^ = V^-. (73) 

For Sk = 1/2, V k these equations were found in [38J through a different technique. 

Through (|72jl we can now determine the energy eigenvalues of ()69j) . It is useful to note 
the following identities 

M M C 

VjS k 



2 Q 5>, + 2££-^ = M(M-l) 
j=i j=i k=i i k 

M C 



3=1 k=l 



Vj - e k 



M C M C 

VjS k v— v v— v s k e k 



EE^-EE^ = ^E- 

j=\ k=l 3 K j=l k=l K k=l 



Employing the above it is deduced that 
c c 

3=1 3=1 
C C C C CM 

E e A- = 2aj2^j + J2J2 s ^- 2M J2 Sk ' 2a Y, v ^ +M ( M -^ 

3=1 3=1 3=1 k^j k=l j=l 

which, combined with the eigenvalues 2sj(sj + 1) for the Casimir invariants Cj, yields 
from (|fj8|72jl the energy eigenvalues 

M c 

E = 2^v j -2^s k e k - (74) 

3=1 k=l 

From the above expression we see that the quasi-particle excitation energies are given 
by twice the Bethe ansatz roots {u,} of ()73|). Finally, let us remark that the eigenstates 
obtained in taking the quasi-classical limit assume the form 

M / C 



e^- io) 



i=i \j=i 

where {vi} satisfy (|73*jl . which are the same as those obtained by Richardson [3H] in the 
case of the reduced BCS model. 

An alternative approach to the exact solution of the reduced BCS model, which pro- 
duces both the eigenstates and the eigenvalues of the conserved operators, was given by 
Sierra |l()()j using conformal field theory given by the SU (2) &-WZW model in the limit 
when the level k approaches —2. There also exists an intriguing analogy for the reduced 
BCS model from two-dimensional electrostatics |101|. I1U2] . 
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7 Exact calculation of form factors 



Through use of the Slavnov formula (|4QI42j) and (j49l50p we have explicit determinant 
representations for the form factors of A(u), B(u), C(u) and D(u). Given any operator, 
we would like to be able to express it solely in terms of these operators, which we call the 
inverse problem. Solution of the inverse problem then permits us to determine the form 
factors for that operator. In the following we will show in several examples how this can 
be achieved. In some cases we will restrict our analysis to some subclass of the models 

(HEED- 



7.1 Form factors for the model of two Josephson coupled Bose- 
Einstein condensates 

Specialising the Slavnov formula ()4UI42j) to the case of the Hamiltonian gives the 
matrix elements of F and G as 

«, - (« - (- + «■» n - n 

In order to apply the Slavnov formula for the computation of wavefunction norms, we 
need to determine the functions (,{v), introduced in (JHJ). We can write 



= C( Vl ,l)---C(v N ,N)\0) 

= (4)^10) 



k=0 

for some scalar functions X{. We deduce from the explicit form of C(u,j) that 

-TV 

X = 7] 

N 

%n = | J (Vi - n.V - (3) . 

i=X 

On the other hand we have 



\v} = B( Vl ,iy---B(v N ,Ny\0) 

N 



k=0 



with 



N 



yo 



Y[ ^ + aN + (3) 



Vn = V 



i=i 

N 
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Using the identity (J35J) applied to the present case 

N 

nK 2 - [«n + pf) = V - 2N 
i=i 

gives 



N 



X N 



V N H(vi-aN-(3)y 

8=1 

N 

V N Y[(v i -aN-P)y N 



which shows that 



i=i 



N 



C(v) = V N l[(v t -aN-f3). 



i=i 



The square of the wavefunction norms are then given by 

Hull = {v\v) 

= C(v)S(v : v) 

where S(v : v) is expressible in terms of F or G. 

From (J2fij) we can see that solution to the inverse problem is achieved through 

a\ = lim — C(u), 

u— >oo U 

a 2 = lim —B(u). 

Using the Slavnov formula we have, for \v) , \w) both eigenstates of the Hamiltonian, 

(v\ ai \w) = (w\ a\ \v) 

= C(w) {w\ a\ \v) 

= ((w) lim — {w\ C(u) \v) 

u— >oo 11 

= ((w) lim -S(w : v ± , ■ ■ ■ ,v M -i,u) 

u— >oo IX 

where 



_Gy = G y for j ^ M, 
GiM = V ■ 



In a similar way we find 



(w \ a 2 \ v) = (v \ a,2 \ w) 

= ((v) {v\ a 2 \w) 
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= £(if) lim — {v\ B{u) \w) 

u—>oo u 

= ((v) lim — {w\ C(u) \v) 

u—>oo U 

= ((v) lim -S(w : v h ■ ■ ■ , v M -i, u) 

with G as above. 

In Josephson's original proposal (HU E] for tunneling of Cooper pairs through an 
insulating barrier, the effect is a manifestation of the relative phase difference of the 
wavefunctions for the two superconductors. Josephson exploited the fact that the BCS 
variational wavefunction ((7|) is not an eigenstate of the total particle number, and as 
phase and particle number are canonically conjugate variables, a well defined relative 
phase could be assigned. For the model (Q) there are technical difficulties which prevent 
a simple definition for the phase variable jlHl [103 • Consequently, the expectation values 
for the Josephson tunneling current 

J = i(a\d2 — c4 a i) 

as well as h and n 2 , where n = Ni — N2 is the relative particle number operator, are of 
primary interest. In principle, these can all be expressed in terms of the form factors for 
a\, a{, a 2 and a\ through completeness relations. This would yield expressions comprised 
of sums of determinants. However, in the case when 

Uu = U 22 , Mi = ^2, (75) 

which results in a = (3 = from (J57|) . we can use a direct method to yield the form 
factors for J~, n and n 2 , expressed as single determinants |H21- The reason we can achieve 
this under the constraint ()75|) is that in this case the Hamiltonian acquires the additional 
symmetry 

[P, H] = 

where P is the permutation operator defined by the action on the Fock basis 

R(a\y(4) k \0} = (a\) k (aly\0). 

This means that the energy eigenstates are also eigenstates of P, and moreover, P 2 = I 
shows that P has eigenvalues ±1. Only by exploiting this symmetry do the form factors 
for h, h 2 and J become accessible. 

As mentioned earlier, the realisation of y[g/(2)] used to derive the model (Q) is not 
unitary. It is however equivalent to a unitary representation when ()75|) is satisfied in the 
sense that 

C\u) = PB(u*)P. 
Consider for {v{\ satisfying the Bethe ansatz equations 

(v\v) = ((v)S(v : v) 

= ((v) (0\B(v N ) ■ ■ ■ B{ Vl )C{ Vl ) ■ ■ ■ C(v N )\0) 

= ((v) (0\PC\v* N )P..PC\v* 1 )PC(v 1 )...C(v N )\0) 

= ((v) (0\C\v* N )...C\v* 1 )PC(v 1 )...C(v N )\0) 

= ((v)(v\P\v) 
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which shows that ((v) = ±1 is the eigenvalue of P; viz. 

P \v) = C(v) \v) . 

From the Slavnov formula, the squares of the norms of the eigenstates in this limit 

||#|| 2 = (v\v) 

= \S(v:v)\ 

are obtained directly. 
We define 

S = A(0) - D(0) = rfNtNz + V~ 2 - 

Letting \v) and \w) be eigenstates of the Hamiltonian we can appeal to and (|S0|) to 
find 

{am = -c^f'^,^ det + 2Q) (76) 

where the elements of Q read 

T 3 Yl?=i( v j ~ v i -v) 



Qi 



Wi(Wi + 77) 



We remark that because the basis states are also Hamiltonian eigenstates, it is straight- 
forward to write down the time-dependent form factors 

(w\E(t)\v) = exp{-it {E{w) - E(v))) (w\E\v) (77) 

where the energies are given by ([58)1. with a — (3 — 0. 

Remarkably, from equation (fTUj) all the form factors for h, h 2 and J can be obtained. 
This is achieved by exploiting the symmetry of the Hamiltonian under P. We begin with 
the following result, which is easily proved. If C(^) 7^ CC^O then 

(wlNiNilv) = 0. 

If ((w) = ((v) then 

(w\J\v} = 0. 
The result follows from the observation 

PN X N 2 = A^iiV 2 P, 
PJ = -jp. 

We now find that 

(w\NiN2\v) = rf 2 {w\B\v) + i]~ 4 (w\v) 
if C(^) = C(v)i an d is zero otherwise. Also 

(w\J'\v) = —i (w\E\v) 
if C(^) C(v)> an d i s zero otherwise. 
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The above shows that the form factors for J are obtained directly from those of 5. 
Those for n 2 also follow, since we have n 2 = N 2 —4N1N2 and the Hamiltonian eigenstates 
are also eigenstates of the number operator N. Thus 

(w\h 2 \v) = N 2 (w\v) — 4 (w\N 1 N 2 \v) 

= (N 2 - 4r]~ 4 ) (w\v) - 4r]- 2 (w\Z\v) 

if C(w) = ((v), and zero otherwise. To obtain the form factors for h, we use the fact that 
J is the time derivative of n, so 

J=i-[h, H] 
w 



which gives 



(w\n\v) = (w\J\v) 

\w\z\v) 



if C(^) 7^ C(^) anc i zero otherwise. 
The expectation values 



E w) - E v 



w. < * |9| * > 



(*l*> 

where 9 = h, h 2 or J , and |^) is an arbitrary state, can be expressed in terms of the 
form factors through completeness relations, in a time dependent fashion. In particular, 
for a given \^f) the quantum fluctuations of the relative number operator 

A^;n) = (n 2 )^(n)l 

can be computed from these results. 

The extension of these results to the general case without the imposition of the con- 
straint (J75|) remains an open problem. 



7.2 Form factors for the models of atomic-molecular Bose-Einstein 
condensates 

As both models for atomic-molecular Bose-Einstein condensates are derived from the 
same L-operator (|27jh we may treat the two models simultaneously. In analogy with the 
previous model, we can deduce 




c' = i] lim —C(u) 

it — >oo U 

K + = -r]- 2 lim - (B(uY + C(u)) . 

u — >oo U 

This leads to the following form factors when {wj} and {vj} both satisfy the Bethe ansatz 
equations 

(v, k\c\w, k) = (w, k|c' \v, k) 
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= r] 1 C(w) lim -S(w: v a , • • • , f M _i, u), 

v, — >oo u 

(v, k\K~\w, k) = (w, k\K + \v, k) 

= -V~ 2 (C(w) + C(#)) lim S(w : t>i, • • • ,v M -i,u). (78) 

u— >oo 

Note that for this class of models we include the label k in the Bethe states in order to 
identify the pseudovacuum used for the Bethe ansatz calculation. 

Realising the su(l, 1) algebra in terms of the Heisenberg algebra as in (|60|) or (pj3"j) 
gives form factors for the models (JHJ) and (J1J respectively. Certain form factors for the 
single particle atomic creation and annihilation operators can also be obtained by using 
the fact that for these models there are multiple possible pseudovacuum states for the 
Bethe ansatz calculations. For the model © we have 

(v,l/4\a\w,3/4) = (w, 3/4\<J\v,l/4) 

= <«?,3/4|iT,3/4>, 
(v,3/4|a|w,l/4) = <w,l/4|aV,3/4) 



(w, 1/41^1^1/4) 

2(td,l/A\K + \v,l/A) . (79) 



'2k 
1 



In the case of the model (@J) we find 

(v, (k — l/2)\a\w, k) = K,\a* \v, {k — 1/2)) 

= \/2k — 1 (w, k\v, k) , 

(v,(K+l/2)\b\v},K) = (w,K\tf\v,(K + l/2)) 

1 

[w, K\a 



w,k\K + \v,k). (80) 
V2k 1 

Note that in the case of (f?9*jl . if |iT 3 3/4) is an eigenvector of ©, there is no reason 
to assume that \v, 1/4) is also an eigenvector. Hence the formula (fTHj) cannot be used to 
evaluate (fT9|) . since in (|T8|) it is required that \v, 1/4) is an eigenvector. This is because 
we can only establish that (pH^ holds for eigenvectors. A similar situation applies to (jEDJ)- 

In the quasi-classical limit the procedure for computing the form acquires a simplified 
form, in which (|79I80|) can be evaluated. Moreover, the form factors for K z can be 
obtained, which are seemingly intractable in the general case. Below is a detailed account. 
We set the coupling parameters Uij to zero in the Hamiltonians (13|4[) . This corresponds 
to the ideal gas limit in the sense that the terms with coupling Uij describe the S'-wave 
scatterings between the particles. Mathematically, this means that r] = 0, corresponding 
to the quasi-classical limit, and uj(x) = (3 is constant. We scale the generating elements 
A(u), B(u), C(u), D{u) of (|27)l by a factor of 1/u, and in taking the quasi-classical limit 
we obtain the following realisation of the Gaudin algebra 

A{u) = —-(u + /3)I, 
u 



K~ 



B{u) 

u 

C(u) = c f 

u 

K z 



V(u) 



u 



39 



with ^ ^ 

a(u) = U — (3, d(u) = . 

U U 

This realisation is evidently not unitary, but it is clear that 

(w, K \ = (-1) M {w,k\ 

even for arbitrary {wi}. Through using the quasi-classical limit of the Slavnov formula 
(15 1|) we may find the scalar product of the states in these models 



(w, k\v, k) = (—1) {w, k\v, k) 

1 



detjF 



Uk>l( v k - vi)Ui<j( w i - Wj) 
M n - n 'f£r' t, ' > r detg (82) 



where we have defined 



T = —T 

% 3 l 3 



alVn) — (C[Va 



2 \ Lh=i( v 3- W i) 



u 3 U -\ U 3J ' / j I I \o 

2 \ '0 ' "V) 



Vj-WkJ (v J -w l ) 2 

in order to absorb the factor (— 1) M . The set {w^ provide a solution to the quasi-classical 
limit of the Bethe ansatz equations 

2k 



Z,K v— v ^ 

Wi + 13 = V . 83 

In the quasi-classical limit the energy eigenvalues for the Hamiltonian (jSJ) are 

M 

E{v) = Ha(2M + 2k - 1/2) - 2Q ^ ■ 
while for (J1J) they are given by 

E(v) = i^ a (M + 2k - 1) + /i b M - Q 



M 

U 



M 

'w, 

i=l 



In deriving the above energy expressions we have used the identity 

M M 



1=1 1=1 



which follows from (jHSJ)- 
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For the case when {vj} also satisfy the Bethe ansatz equations we find for the elements 
of T 



Vj -W k I (Vj - Wi) 



2 



(I l\ ^ 2 ™ 2 
[2k( + Wi - vj + V V 



k^ti J k=f=i 

x uZi( v j- w i) 

2k M 2 \ M 



W i V 3 TZl ( W i ~ W k)( v j - W k ) 



and similarly 



k^i y 1 "7 l^i 



e, = (i + — + ti 1 rlr^- < 85 > 

\ T+i ( Wi ~ W k)( V 3 - W k) J ( v j - w i) 

Letting \v,k) = \w, k) gives us the square of the norm formula 

I \v\ I 2 = det /C 



where 



—7 2k ^ — > 2 



2 

(U< - Vj)' 



fcij = 'Jr. ITT? for * ^ 



To compute the form factors for K z , we need to take a limit of the form factors of 
T>(u) as given by This leads us to 

(u>, k) = — lim u (w, k\T>(u)\v, k) 

«— >0 

T-rM 

K \h=l W k 



m i m M, m M det(^-Q) 



Yi k >i( v k - vi) - ^i) ru=i 

with 

~ „.,2 



"7 



Using the fact that K z + N c is conserved in both models ()3I4|) and can be expressed in 
terms of the total atom number N, the form factors for N c can be deduced from those for 
K z . 

Next we turn to the problem of finding the form factors for the operators c\ c, K + and 
K~ . To do this we need to solve the inverse problem and express each of these operators 
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in terms of the realisation of the Gaudin algebra. This is not difficult to achieve with the 
result 

K + = — lim uC{u), 
K~ = lim uB(u), 

u— >0 

= lim C(u), 
c = - lim B{u). (86) 

u— »oo 

Using the fact that the parameters {vi} in the Slavnov formula (|82|) are arbitrary, we 
can then take the limits described above to yield the form factors. This gives the results 

(v, k\K~\w, k) = (w,k\K + \v,k) 

= — lim u (w, k\vi, ■ ■ ■ , vm-i, u, k) 

rM 



~\~\1V1 

U '=> W < ietV 



where 



u k >i (vk-vi)Ui<j\«>i-wj)np=i v p 



Vij = T %3 for j/M, 
2k 

2 ' 



and 



(v,k\c\w,k) = (w,k\c*\v,k) 

= lim (w, k\vi, ■ ■ ■ , v M -u u , K ) 

n - iuk-".) detw 



where 



Wij = g tj for j^M, 
W tM = 1. 



7.3 Form factors for the reduced BCS model 

The results of this section have been published in [HI] for the case Sk = 1/2, VA; (although 
different conventions and notations were used). A closely related study is given in |l()4j . 
The fundamental difference between jlU4j and the results below is that |104j employs 
the generating function of correlators of the Gaudin algebra as developed in [3T] , whereas 
below we will directly use the quasi-classical limit of the Slavnov formula as given by (|5ip. 
By this procedure the form factors are obtained in an explicit determinant representation. 
Again, we will derive results for the general case of the Hamiltonian where the 
irreducible realisations of the su(2) algebras, labelled by a lowest weight — Sk, are arbitrary. 
The realisation of the Gaudin algebra obtained by taking the quasi-classical limit of the 
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realisation of the Yang-Baxter algebra given by reads 



A{u) = + ' 



B{u) = Y-^. 

an) = J2— 



c Si 



k=l 

with 



V(u) = aJ-y;— (87) 
^ u - e k 



c 

a(u) = — a — , d(u) = a + ^^ 



, -. u - e k u € k 

k=l k=l 

This realisation is unitary, so we do not need to deal with the issues of non-unitarity as 
in the previous examples. 

Through using the quasi-classical limit of the Slavnov formula (|51|) we find the scalar 
product of the states 

(w\v) = — jj jj detJF 

Uk>l( V k ~ Vl) UiKj K - Wj) 



nM t-tM , \ 
p=ill 9 #pK-^g) 

Uk>i( v k - vi) ni^-K - w i) 



det£ 



where Qij and T%j are given by (|52I53|) and the Bethe ansatz equations for the parameters 
{wi} are given by (f75|). Letting {v^ also be a solution of the Bethe ansatz equations we 
find 



2sfe 2 ^ (vj - wj) 

• • in- — £i 



fc= ^ (wi - e k )(vj - e fc ) ^ - w fc )(u,- - w fe ) y (fj - ^) ' 

Setting \v) = \w) gives us the square of the norm formula 

I \v\ I 2 = det K, 



where 



k=l fc^i v 

fCij = — for i ^ j. 
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Putting Sk = 1/2, V/c, the above is exactly the norm square formula obtained by Richard- 
son [THE] , 

To compute the form factors for S^, we use the fact that T>(u) has simple poles at 
u = ej, Vj; i.e. 

S z m = - lim (u - e m )V(u). 

This leads to 

(w\S^\v) = - lim (u - e m ) (w\V{u)\v) 



where is as above and 



det(JF - Q( €m )) 



Qn(u 



%3 (u - Wi) 2 

Now we derive the form factors for the operators S + and S~. In this instance the 
inverse problem is solved as follows 

S~ = lim (u- e m )B{u), 

S+ = lim (u-e m )C(u). (89) 

Using the fact that the parameters {v i} in the Slavnov formula (jHHj) are arbitrary, we can 
then take the limits described above to yield the form factors. The results are 

(v\S~\w) = (w\S+\v) 

= ]im(u-e m ){w\v 1 ,---,v M - 1 ,u) 



nd X (% - e m) Uk>l\ V k - Vl) Ut<j( W i ~ W 3 



detP 



where 



Pij = Fij for j^M, 

V iM -- ' 



(■Wi - e m ) 2 



The above form factors can be used to construct general correlation functions, such 
as the Penrose-Onsager-Yang off-diagonal long-range order parameter as given in [oTj . 

8 Conclusion 

We have reviewed the theory of the quantum inverse scattering method and algebraic 
Bethe ansatz for the computation of energy spectra and form factors in exactly solvable 
models, and demonstrated how it applies to several models of Bose-Einstein conden- 
sates and the reduced BCS model. Throughout we have only used the specific example 
of the Yang-Baxter algebra associated with the Lie algebra gl{2). However, a Yang- 
Baxter algebra can be associated with any simple Lie algebra, Lie superalgebra, and the 
g-deformations of these structures. Hence the theory can be applied on a much wider 
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level. For example, generalised BCS systems derived from Yang-Baxter algebras associ- 
ated with the Lie algebras gl(4) and so(5) and the Lie superalgebra gl(2\l) were derived 
in |lUb't HU7t HU8j respectively. The model obtained in |lU7j was proposed by Richardson 
in 1966 to describe proton-neutron pairing in nuclear systems |109t 1110"] . The case of a 
general Lie algebra was examined in in the context of the Knizhnik-Zamolodchikov 
equation. One major challenge which remains is to extend the Slavnov formula for the 
scalar products of states to the general case. In fact there has been little progress on this 
aspect with the exceptions of the work by Reshetikhin on the norms of the wavefunctions 
for models derived from the gl(3) Yang-Baxter algebra |112j and by Gohmann and Korepin 
on the Hubbard model |113| . Another approach based on the Knizhnik-Zamolodchikov 
equation for gl(n) can be found in j!14j . 

In the method we have described the exact solution is parameterised in terms of the 
Bethe ansatz equations, which cannot be solved analytically. Consequently numerical 
analysis of the solutions must be undertaken. For the BCS model there have been quite 
a number of works on this topic (e.g. see (10J I102[ ll(J4[ I115| I116j ). For the models of 
Bose-Einstein condensates the only numerical analysis of an exact solution, of which we 
know, is in jnHj. It is also possible to conduct an asymptotic analysis of the Bethe ansatz 
equations to find the exact asymptotic behaviour of the energy spectrum and correlation 
functions. Examples are given in JHH1 EH EH3 ■ 
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